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ABSTRACT 

We argue that the first stars may have spanned the conventional mass range rather than be iden- 
tified with the Very Massive Objects 100 — IO^Mq) favoured by numerical simulations. 
Specifically, we find that magnetic field generation processes acting in the first protostellar 
systems suffice to produce fields that exceed the threshold for MRI instability to operate and 
thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on proto- 
stellar mass scales below ~ SOM©. Such fields allow primordial star formation to occur at 
essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion 
and feedback in much the same way as occurs in conventional molecular clouds. 

Key words: stars: formation: general - stars: primordial - galaxies: star formation - cosmol- 
ogy: magnetic fields - metalHcity 



1 INTRODUCTION 

The formation of the first stars should in principle be simpler to 
understand than present-day star formation. The usual prediction 
that the mass of the first stars 0(100 - 1OOO)M0 is attributable 
to dominance of H2 cooling and consequent high temperature and 
accretion rate (Abel, Bryan and Norman 2002; Bromm, Coppi and 
Larson 2002). It is appropriate to question this result, that only very 
massive stars formed, for at least four phenomenological reasons. 

Firstly, one finds solar mass stars at [Fe] — 4. These may 
be contaminated by companion stars that exploded as hypemovae, 
whose ejecta are depleted in Fe relative to N and C and this in- 
terpretation is consistent with the measured abundance ratios for a 
handful of extreme metal-poor stars. However not all metal poor 
halo stars display such anomalies, and a substantial fraction at the 
10~* solar level have abundance ratios that are consistent with con- 
ventional SNII precursors. There are at least two halo stars known 
at [Fe] < - 5. HE 0107-5240 is a giant with [Fe/H] =-5.3, but en- 
hanced nitrogen, carbon and oxygen: [N/H]= -3.0, [C/H]= -1.3 and 
[0/H]=-2.9 (Christlieb et al. 2004). HE 1327-2326 is a main se- 
quence star (or subgiant) with an iron abundance about a factor of 

2 lower than HE 0107-5240. In this latter case, both nitrogen and 
carbon are enhanced relative to iron by about 4 dex, while there is 
only a comparable upper limit on oxygen (Frebel et al. 2005). 

Appeal to preenrichment by a core collapse ~ 25 M© super- 
nova of Population III abundance with fallback fits the abundance 
patterns well (Umeda and Nomoto 2003), although oxygen may 
possibly be overproduced in the case of HE 0107-5240 (Bessell, 
Christlieb and Gustafsson 2004). The example of HE 1327-2326 
eliminates any internal enrichment source (e.g. convective dredge- 
up), nor is there evidence for a binary companion in the case of HE 
0107-524, thereby also suggesting that mass transfer from an AGB 
star is an unlikely explanation. The high CNO abundances argue 
for a common explanation involving enrichment of the primordial 



cloud by Type II supernovae of primordial abundance. Indeed, the 
abundance patterns in extremely metal poor halo stars suggest that 
enrichment was produced by Population III stars in the mass range 
20-130 A/0 (Umeda and Nomoto 2005). 

Secondly, the broad emission line regions of very high redshift 
quasars, reveal high elemental abundances that appear to have also 
been generated by conventional SNII precursors. In particular, the 
nearly constant Fell/Mgll emission line ratios over < 2 < 5 
requires intense SNII activity at a redshift as high as 9 (Dietrich, 
Hamann, Appenzeller and Vestergaard 2003). 

Thirdly, the chemical yield predictions from primordial very 
massive stars, when normalised to the inferred ionising photon out- 
put required to reionise the universe at high redshift, do not cor- 
respond to observed abundances in any primitive environments. In 
particular, the pair-instability supernova nucleosynthetic signatures 
generated by stars with initial masses in the range 130 to 260 Mq 
are not seen either in the intergalactic medium, including both Ly- 
man alpha forest and damped Lyman alpha absorption systems, or 
in extremely metal-poor halo stars (Daigne et al. 2005). 

Finally, from the theoretical perspective with regard to cool- 
ing, the primordial metal abundance pattern has profound conse- 
quences for the thermal balance and chemical composition of the 
gas, and hence for the equation of state of the parental cloud. 
Spaans and Silk (2005) find that the polytropic index is soft for 
low oxygen abundance enhancements, [0/H]< —3, as appropri- 
ate for Population III, but stiffens to a polytropic index 7 larger 
than unity for [0/H]> 10"'^ due to the large opacity in the CO 
and H2O cooling lines. Hence Pop III star formation is efficient, 
especially before hypernova enrichment and associated oxygen en- 
hancement has occurred. There should be no obstacle to forming 
stars over a wide range of masses even in the absence of significant 
fine-structure cooling. On the contrary, once the polytropic index 
stiffens, at [0/H] — 2, the IMF should change and, specifically. 
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flatten, as argued by Spaans and Silk (2000). Star formation is sub- 
sequently less efficient, in part because of massive star feedback in 
addition to the stiffening of the equation of state, as is required in 
most discussions of galaxy formation in order to avoid premature 
consumption of the gas reservoir. 

Of course, to do justice to the very massive primordial star 
h5rpothesis, the typical expectation of one very massive star per pri- 
mordial cloud is unlikely to produce enough metals to raise the 
mean IGM metallicity to a value of ~ 10~* — 10~^'^ solar, as ob- 
served in the Lyman alpha forest (Schaye et al. 2003) and as also 
required to significantly modify the cooling and thereby enhance 
fragmentation. This is especially plausible as the ejecta are mostly 
not retained in the shallow potential wells of the pregalactic clouds 
(Norman, O'Shea and Paschos 2004). Rather, one has the impres- 
sion that the first stars that produced metals in significant amounts, 
in particular as monitored in the abundance patterns of the Lyman 
alpha forest at modest overdensity and in halo stars at [Fe] ^ — 4, 
spanned the conventional stellar mass range. 

A more fundamental issue may be that the elegant state of the 
art simulations with cosmological initial conditions have nothing 
directly to say about the primordial stellar initial mass function, nor 
for that matter about the physics underpinning the inferred trans- 
fer of angular momentum, but only about the mass function of gas 
clumps. The ultimate fate of these gas clumps has not yet been con- 
vincingly demonstrated. 

Notice that recent simulations by Saigo, Matsumoto and 
Umemura (2004) show that Population III stars are prone to form 
in binary systems through the formation of a rotationally supported 
disk. For that to happen, the initial angular momentum of the par- 
ent cloud needs to be small, the centrifugal force being at most 
10% — 30% of the pressure force. However, the accretion efficiency 
seems not affected by the binary formation, and the binary Pop III 
stars are again predicted to be very massive. 

Of course, sdl of this begs the question as to whether a soft- 
ened EOS can lead to low or intermediate mass star formation in 
situations of very low metallicity, or at any rate, reduce the charac- 
teristic stellar mass to below the lower limit for pair instability SN. 
Certainly, the clouds from which stars with [Fe/H] — 5 formed 
seem to have had a soft equation of state. The minimum fragment 
mass is lowered at [O/H] ^ — 3 due to enhanced fine-structure 
cooling (Bromm and Loeb 2003), but it is likely that turbulent frag- 
mentation (Klessen et al. 2004) enhanced by the softening of the 
equation of state is the key to low-mass fragmentation, for example 
via disk formation. According to this viewpoint, it is the equation 
of state in the initial collapsing cloud that controls the IMF. 

An alternative view is that the final fragment mass remains 
high, as suggested by the numerical simulations, but that a second 
mode of fragmentation associated with the onset of dust cooling at 
very high density may lead to low mass star formation (Omukai et 
al. 2005). The latter study finds that a metallicity of Z ^ IQ^-'Zq 
at a density in excess of 10^"cin^'^ marks this transition. How- 
ever Spaans and Silk (2005) argue that the gas properties, and in 
particular the equation of state, at much lower densities, typically 
the range associated with saturation of H2 cooling at ~ 10*cm~^, 
most likely determine the IMF. 

The purpose of this paper is to demonstrate that a plausible al- 
ternative to dust cooling in primordial clouds involves recourse to a 
mechanism that is universally accepted as being an essential key to 
controlling the IMF in conventional molecular clouds at the onset 
of collapse, namely magnetically-regulated fragmentation and ac- 
cretion. The essence of our argument will be to derive the magnetic 
field strength in primordial star formation. An important conclu- 



sion is that one should continue to find low mass, and especially 
abundance signatures of intermediate mass, stars down to essen- 
tially zero metallicity. Note that the possible effects of magnetic 
fields on primordial star formation have also been investigated in a 
pioneering paper by Tan and Blackman (2004), and more recently 
by Machida et al. (2006). 



2 THE TRANSITION FROM PRIMORDIAL STAR 
FORMATION 

For the moment, pending higher resolution simulations with ade- 
quate exploration of the EOS, it may be grasping at straws to argue 
that Pop in could contain a substantial or even a dominant low mass 
component. Nevertheless, an analytic exploration is warranted. In 
this note, therefore, an alternative approach is addressed. A reason- 
able hypothesis, prevalent in the star formation community, is that 
magnetic fields are essential to the generation of the initial stel- 
lar mass function in conventional star-forming clouds (e.g. Mestel, 
1965; Lizano and Shu, 1989; Shu et al., 2004), even if precise de- 
tails of the actual regulation mechanism still remain to be worked 
out. If this is indeed the case, one may then ask what happened 
in the first clouds, when presumably magnetic fields were absent 
at any significant dynamical level? We refer to this as Primordial 
Star Formation. And was there a seamless transition to Current 
Epoch Star Formation, involving no drastic change in the under- 
lying physics? 

One common, and generally accepted, viewpoint among cos- 
mologists is that metallicity was a primary driver When the metal- 
Ucity was less than about 0.001 of the solar value, the dominant 
coolants were H atoms and H2 molecules. These sufficed only 
to cool to about 200K. The inefficient cooling is the primary rea- 
son that both analytic and numerical calculations of the minimum 
fragmentation mass in primordial clouds yields high accretion rates 
onto protostellar cores (M* w /G) that result in characteristic 
primordial stellar masses of ~ 100 — lOOOM© . 

However the protostellar core is only O.OIMq and grows by 
accretion. One may try to understand how low mass stars could 
form in a primordial, or nearly primordial, environment if frag- 
mentation, itself inefficient, was aided by magnetic feedback. Con- 
sider the following hjfpothesis. Magnetic fields are generated by 
a dynamo associated with magneto-rotational instabilities at work 
in protostellar or circumstellar disks. There are several necessary 
conditions for such an MRI dynamo to operate. Firstly, the proton 
Larmor radius must be less than the disk scale-height. Secondly, 
there are important conditions imposed on the wavelength and on 
the growth rate of the MRI instability. Thirdly, there must be a seed 
field that exceeds the minimal field under which the MRI d3mamo 
operates. In the following, we will examine each one of these is- 
sues, within the framework of a simple steady state circumstellar 
disk model. This will allow us to obtain a critical mass of the cen- 
tral (proto-)star above which all the conditions are fuUfiUed for the 
MRI to work. 

2.1 General equations 

We start by considering a three fluid model, including free elec- 
trons, ions (protons, for simplicity) and neutral gas, of respective 
number densities 11^,111 and rin- Setting n = rie = n-i global 
charge neutrality implies that the electric current is j = —enV, 
where V = Ve — Vi is the electron velocity relative to ions. Simi- 
larly to the procedure followed by Mestel (2003), we can combine 
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the equation of motion of electrons witli tliat of ions to obtain tlie 
following generalized Ohm's law equation, 
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where k~ — lot, — ujTe and — uiTi measure how 
strongly the charged particles are coupled to the magnetic field 
with respect to their collisional interactions (oj = eB/m^c and 
uji = eB jvnic are the electron and ion gyrofrequencies respec- 
tively, r being the mean electron-ion collision time, and and 
Ti the mean times for electron and ion collisions with neutrals). 
Note that the products KqB are independent of the magnetic field 
amplitude B. The density fraction of neutral particles is given by 
F = pn/ [pi + pn), and their fractional contribution to the total 
pressure p is / = Pn/p- The rest of the symbols used above bear 
their usual meanings. 

We get the magnetic field induction equation by taking the curl 
of Eq. G}. I gnoring the pressure terms, we obtain 



dtB = V X 




(2) 



in a local Keplerian frame corotating with the disk at angular fre- 
quency Q,. The second, third and fourth terms on the right hand-side 
are the resistive. Hall and ambipolar diffusion terms respectively. 
It is straightforward to check that we recover the usual induction 
equation in the limit of vanishing density of neutrals (F, f — > 0, 
and Ke,i since Te,i — » oo). 

Similarly, the equation of motion of the whole plasma is 

= pV(t> + + F,, - p{2Q.x V + Q.X {q.x?^^ (3) 

where F^ accounts for the viscous forces, and d- /dt — df+Clx- + 
i7- V- is the full time derivative. Let us define the velocity departure 
from Keplerian motion hy Sv = v — Q,reg. The equation of motion 
then reads 

+ ^dg + 5iJ ■ vj (Su = —2Q,{5vreg — Svge,.) + Sfi^reV 
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and the induction equation simply becomes 

dtB = V X '^{Sv + Q.reg) x B 

-V X \ Ke + K, — ) J 



B 



(5) 



where we omitted the Hall and the ambipolar terms which are of 
second order for initially weak magnetic fields. 

Circumstellar disks are prone to gravitational instabilty. Es- 



pecially in the context of Pop III star formation. Tan and Black- 
man pointed out that proto-stellar disks may develop gravitation- 
ally driven turbulence, amenable to the a formalism (e.g. Gammie, 
2001), that provides angular momentum transport. As shown by 
Eqs. J4j and J5}, gravitational instability may amplify velocity fluc- 
tuations, which in turn can, in principle, drive the growth of mag- 
netic fields. The importance of gravito-turbulent effects depends on 
the ratio of the first term and the last term of the right hand side 
of Eq. j5j, which yields the actual magnetic Reynolds number. The 
factor multiplying the current on the right hand side of Eq. J5} is 
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Taking the mean collision time of charged particles with neutrals 
from Draine et al. (1983), we obtain 



Te \ 104K 
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and 



Te 



and the magnetic Reynolds number is then 
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where the square brackets stand for the last term in equation ^6), 
accounting for the effects of neutrals. The presence of neutral par- 
ticles in the disk would enhance the dissipation effects. For initially 
weak magnetic fields, this would reduce the capability of gravita- 
tional instability to amplify any weak magnetic seed fields, even if 
velocity fluctuations would be amplified. Taking the velocity fluc- 
tuations to be of the order of the sound speed in the disk, and the 
gradient scale being at most of the order the disk radius, the mag- 
netic Reynolds number reaches ~ 10^* at T ~ 10^ K. 

Based on the high Reynolds number. Tan and Blackman 
(2004) argue that local gravitational instability develops, leading to 
an alpha-driven dynamo. This could provide a mean of amplifying 
any initial seed fields to have a dynamical effect on primordial star 
formation. However, recent simulations of massive self-gravitating 
disks find that global gravitational instability modes dominate the 
energy transport (Lodato and Rice 2005). If this is correct, as Bal- 
bus and Papaloizou (1999) initially suggested, the alpha formalism 
would not be effective in field amplification. 

Amplification of magnetic fields may however be driven by 
the action of MRI turbulence. The role of neutrals, Hall effect and 
ambipolar diffusion in weakly ionized disks has been investigated 
in several papers (e.g. Wardle 1999, Salmeron and Wardle 2003, 
Kunz and Balbus 2004). However, as we discuss below, one does 
not expect proto-stellar disks around forming Pop III stars to be 
cool enough so as to allow a dominant amount of neutral particles. 
In the following, we therefore consider a two fluid disk model only. 



2.2 Simple disk model 

To demonstrate the potential role of magnetic fields in Pop III ac- 
cretion, we use a very simple accretion disk model. We consider 
a thin, nearly Keplerian disk in the steady state. In that case, the 
accretion rate M is constant throughout the disk, and the surface 
density 
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can be related to the accretion rate via 
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wiiere p(po) is the disk (midplane) gas density, H is the disk height 
scale, V is the viscosity, n is the disk inner radius and /3 is a param- 
eter related to the angular momentum flux (see, e.g.. Spruit 2001). 
We rewrite the viscosity in terms of the dimensionless a viscosity 
(Shakura and Sunyaev, 1973) as i/ = acsH, where 



1/2 



(12) 



is the sound speed (jj, ~ 0.6 is the mean atomic weight for a gas of 
primordial composition, and / = H/r is the disk aspect ratio). We 
have used here the relation between the actual temperature in the 
disk and the virial temperature, 

T = /r...l.4xl0V^^(^)"K. (13) 

Combining eqs. (To}, O and (T2j, we obtain an expression for the 
midplane density 
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Furthermore, if the central star rotates at a rate fi* < f^Kcp, one 
finds that (3—1, and at large radii r as compared to the disk inner 
radius r^, the last term between the square brackets above can be 
ignored. 

We are now left with the task of evaluating the aspect ratio, /. 
The few available models of disks around a forming Pop III star do 
not really agree on the disk thickness (and overall structure). Tan 
and McKee (2004) found a thin disk solution with / ^ 0.13 up to 
r ~ 20 AU whereas Mayer and Duschl (2005) show that, in their 
"fiducial model" with M = lO'^Mo/yr, / ~ 0.4 up to r ~ 46 
AU, and grows to reach unity at larger radii. 



2.3 Magnetic field seeds in the disls 

Various models have been proposed for the creation of cosmo- 
logical magnetic fields in astrophysical plasmas in the early ages 
of the post-recombination Universe. Most of them rely either on 
thermal pressure effects (the so-called Biermarm battery, Biermann 
1950; for applications in cosmology, see e.g. Kulsrud et al. 1997, 
Hanayama et al. 2005) or on radiation drag (e.g. Harrison 1970, 
Langer, Puget and Aghanim 2003, Langer, Aghanim and Puget 
2005) operating on cosmological distances. Those usually lead to 
relatively weak fields on large scales, and one has then to rely on 
amplification by contraction or dynamo and make additional as- 
sumptions (e.g. flux freezing, geometry of collapse, etc.) in order 
to get interesting values of fields in circumstellar disks. A particu- 
larly effective turbulent amplification mechanism that operates in a 
dilute plasma such as the intracluster medium or galactic halo ap- 
peals to a self-accelerating fluctuation dynamo where anisotropic 
pressure gradients generate plasma instabilities that result in a high 
Reynolds number (?). However, the same physical processes are 
also at work in the disk itself, and comparatively high fields may be 
directly created in situ. We therefore do not need to have recourse 
to preliminary magnetic seeds generation on larger scales. 



For a two component system, and with rrie <^ pmp, combin- 
ing the ion and electron equations of motion leads to the general- 
ized Ohm's law 
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the curl of which yields the induction equation 
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The first term on the right hand side above corresponds to the dy- 
namo amplification, and since we are interested here in the creation 
of magnetic fields starting with B = 0, we will ignore this term 
in the following. The second and last are the battery and radiation 
drag terms respectively. Both can serve as sources to create a seed 
of magnetic fields in the disk. 

As is well known, the Biermann battery term is non zero only 
when the temperature and electron density gradients are not paral- 
lel, according to 
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The characteristic time scale is the rotation time f2, and the 
electron density and temperature vary over scales at least of or- 
der H and r. Therefore, we expect the field created by the battery 
mechanism to be of order 

ck 1 T 

-Dbatt ~ ,/ — , 
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1/2 
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Turbulence might actually create density gradients on smaller 
scales, leading to higher field amplitudes, but by considering longer 
length scales, we adopt the more conservative approach and keep 
the lower estimate for B. 

The value we obtained from the Biermann effect is just a factor 
/ smaller than that provided by radiation drag. Indeed, as has been 
shown by Balbus (1993; see also Chuzhoy 2004) through momen- 
tum conservation arguments, the typical amplitude we may expect 
is 
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It so happens that this amplitude. Brad, corresponds to values of the 
ions Larmor radius of order the local disk radius. This means that 
the magnetic field is only weakly coupled to the ions. Therefore we 
expect the advection of field lines by small scale turbulence to be 
quite inefficient. 

2.4 MRI dynamo instability 

We now wish to derive a constraint on the mass of the central object 
by requiring the magnetic field in the disk to exceed the minimal 
value for the magneto-rotational instability (MRI) to operate and 
drive a dynamo. 

The first basic condition we consider is that the wavelength of 
unstable modes must exceed the particle mean free path. This con- 
dition gives the minimum field for MRI. Now in MRI (Balbus and 
Hawley 1998), one expects that the dominant growth mode satisfies 
Ami ~ va/^, where va is the Alfven frequency and f2 is the disk 
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rotation rate. Larger scales become unstable as the dynamo oper- 
ates, and the maximum scale is given by the disk scale height. This 
gives the maximum field strength where the dynamo saturates. For 
a sustained field, that does not reverse sign every rotation period, 
one requires a stratified disk, which allows helical turbulence to be 
generated. Helical turbulence results in mean field amplification via 
the MRI dynamo and it is possible that the resulting fields can at- 
tain the equipartition value (Brandenburg 2004). Now the mean free 
path in an ionised gas is 9r^fc^/(7rnee* In A) (e.g., Lang 1999), 
and the minimum field condition is thus 

9 (20) 



n ^ InA Trriee*' 
or 



B>B.Mp = 1.75xlO-«a^/V"^'T^f 



11/4 



M 



In A VIOM0, 
1/2 / „ \ -11/4 



(21) 



where we assumed = fiKcp- Notice that MRI may occur even in 
coUisionless plasmas, as shown by Quataert, Borland and Hammett 
(2002). Their study, however, does not account for possible effects 
due to a finite Larmor radius, which we presumably cannot ignore 
in our case where initially weak fields imply very large Larmor 
radii. 

The second condition we may consider is that MRI will be ef- 
ficient in amplifying magnetic fields only if the amplification time 
scale is shorter than the diffusion time scale. For maximum insta- 
bility scales, the growth rate is of order Q., and the condition can be 
rewritten as 



(22) 



where 77 = 10^^(lnA/20)r cm^/s is the magnetic diffusiv- 
ity. Assuming Keplerian rotation, this condition translates to 



Using eqs. jl3> and jl4L we obtain 
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As noticed by Tan and Blackman (2004) who obtained a similar 
result, this value of the amplitude is quite high for initial conditions 
in primordial star formation which starts in a protostellar environ- 
ment not necessarily magnetized beforehand. However, it does not 
imply that MRI could not be at work in a disk surrounding the 
progenitor of a Population III star. Indeed, as already stressed, the 
minimum value obtained above has been derived by taking into 
account effects of magnetic diffusivity on scales of maximal insta- 
bility, which are small with respect to the disk height. This only 
implies that due to diffusivity, MRI for small-wavelength distur- 
bances is de facto inefficient, unless the amplitude of the field is al- 
ready large. Moreover, for standard values of the parameters, Bmi 
is way above Brad- This means that /mi, the scale of maximal in- 
stability, is much smaller than the ion Larmor radius, and therefore, 
as noticed in the previous paragraph, MRI turbulence is unable to 
amplify the fields on such small scales. 

Nevertheless, larger scales of order the disk (or Larmor) ra- 
dius, may still remain unstable, with the growth rate in that case 



being proportional to the field strength itself (e.g. Balbus 1995), 
even if magnetic fields are dissipated on small scales. This brings 
us to the third and strongest condition that must be satisfied. In a 
differentially rotating, turbulent disk, the growth rate for MRI in the 
weak field limit is 7 ~ VAk, where k is the vertical wavenumber, 
where va «C fl/fc and is the local angular velocity. The growth is 
stabilised by magnetic diffusivity, ri, so that the dispersion relation 
is 7 + rjk^ ~ VAk (e.g. Kitchatinov and Riidiger 2004). Therefore, 
for marginal stability. 



VA 



(25) 



with \ < H — fr. This criterion translates readily into a minimum 
magnetic field 

B„„i = 27r=*/'77p'/V"''-"' (26) 



which yields 
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using again eqs. <13> and <14> above. 



2.5 Critical mass of the central object 

We require now that the minimum magnetic field for the MRI to 
work is smaller than the field created in the disk. As we have seen, 
the smallest suitable field amplitude, Bmri, is given by the marginal 
instability criterion applied to scales of order the disk scale. At the 
same time, in situ magnetic field generation is likely to yield am- 
plitudes of order Brad- This amplitude must be higher than Bmri, 
which happens as soon as 

4/9 / j> \ 2/9 



In A 
20 



M 



lO-^Me/yr 



5/9 



(28) 



The thickness of the disk has the most dramatic effect on the value 
of the critical mass, whereas the actual value of the viscous a pa- 
rameter plays a minor role. If the disk is rather thin, with / ~ 0.1 
for instance, then the mass spans a range of rather large values, go- 
ing from 128 to 357A/0. However, in the absence of metals, the 
properties of the primordial gas are likely to prevent the disk from 
efficiently cooling, and we expect / to be closer to unity. In that 
case, taking / ~ 0.4 as obtained by Mayer and Duschl (2005), the 
critical mass is much smaller, comprised between roughly 4.3 and 
27.8Mq for a ranging from 1 to 0.01. 



3 DISCUSSION 

Magnetic fields have probably played a role in Primordial Star For- 
mation, even if starting from a medium free of magnetic fields. As 
we have argued, radiation drag or thermal pressure effects are able 
to generate magnetic seeds in the disk surrounding the central ac- 
creting stellar progenitor Initially, those seeds are weak too much 
for being amplified by small scale turbulence, or even MRI effects 
on scales comparable to the disk scale. However, as matter accre- 
tion proceeds, the mass of the central object grows, and the gravi- 
tational potential it creates deepens, increasing the rotation veloc- 
ity of the disk. Eventually, the rotation is fast enough so that long 
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wavelength MRI modes become unstable, as the minimum mag- 
netic field for MRI becomes smaller than the field generated in the 
disk. Depending on the properties of the disk, this happens once the 
mass of the central object reaches 4 — 28M0. 

Subsequently, MRI dynamo will be at work and will amplify 
the magnetic field. The field can then rapidly reach dynamically 
important values, and magnetically-driven ejection contributes to 
lower the effective accretion efficiency. The actual model of mag- 
netic winds is beyond the scope of this article, but this suggests that 
feedback is likely to require magnetically-driven outflows, which 
could occur during Pop III formation already when, as we argued 
above, the mass of the stellar progenitor is rather small. Note that, 
in their model of protostellar disk in primordial star formation. Tan 
and Blackman (2004) estimated the power of magnetically driven 
outflows. In their study, the outflow feedback effects reduce the star 
formation efficiency once the protostar reaches roughly 100 Mq. 
Incidentally, concurrent conclusions were reached by Machida et 
al. (2006) who simulated the collapse of a magnetized primor- 
dial cloud in rigid rotation, and the subsequent formation of a Pop 
III star, within the ideal MHD approximation. Depending on the 
(high) initial value (Binit > 10~^G for an initial cloud density 
ric ~ lO^cm"^) of the magnetic field, their simulations show that 
magnetically driven jets develop and effectively reduce the accre- 
tion rate. Further numerical studies with higher resolution seem 
necessary to address both the mass ejection rate and the amplifi- 
cation of magnetic fields in proto-stellar disks for initially weaker, 
even vanishing, magnetic seeds. 

We have argued that the interplay between the two modes of 
star formation, primordial, massive and conventional, involving all 
masses, is controlled by B and not by Z. Effective feedback re- 
quires that of order 10 percent of the gas accretion rate is channelled 
into star formation, with an outflow rate that on the average must 
be of the order of the net star formation rate. Globally, for the star 
forming cloud, one expects that MoutSow ~ Af* ~ O.lMaccretion, 
much as is found in nearby cases of star-forming clouds. 

Moreover, feedback is likely to be responsible for the turbu- 
lent support in clouds that lowers the star-formation efficiency and 
helps to generate the conventional IMF. Thus the onset of cloud 
fragmentation, due eventually but we have argued, not exclusively, 
to the enhanced role of cooling, would allow field amplification, an- 
gular momentum transfer, feedback and low mass star formation. 

Alternatively, suppose we accept the hypothesis that the first 
stars were massive objects. To avoid the empirical objections dis- 
cussed above, one would have to argue that merging and coagu- 
lation of gas clumps resulted in formation of predominantly very 
massive objects, of characteristic mass ^ lOOOM : , whose fate is 
to form intermediate mass black holes with relatively low nucle- 
osynthetic yields associated with their collapse. In this case, the 
seed fields may come from jets and outflows driven by spin-up of 
tiu'bulent accretion disk dynamos as a consequence of accretion 
onto these intermediate mass black holes. One can argue that the 
first generation of primordial clouds, which cooled predominantly 
via Lyman alpha emission, preferentially formed IMBHs, since the 
associated minihalos of mass ~ lO^Mo are promising environ- 
ments for forming IMBHs in view of the high core accretion rates 
(Zhao and Silk 2005). Accretion disks around the IMBHs provide 
promising sites for MRI dynamos. 

Alternatively, the jets may generate magnetic fields via the 
Weibel instability, that although small, amounting to ~ mc/rrip of 
the equipartition value, could still be useful as seeds (Wiersma and 
Achterberg 2004). The IMBH outflows allow one to generate larger 
seed fields required for conventional star formation than inferred 



from the Biermann battery. The associated nucleosynthetic impli- 
cations are modest, because IMBHs have essentially zero yield. For 
any reasonable IMF, there are of course associated Population III 
stars. 

Consider for example the inference from primordial star simu- 
lations that one star of mass ~ 10^ Mq forms per pristine pregalac- 
tic cloud of baryonic mass ~ IQ'Mq, a typical galactic precur- 
sor not dissimilar in mass to dwarf satellite galaxies. The resultant 
local eiuichment would amount to of order 10~^ the solar value, 
and would be diluted by a further factor that corresponded to the 
fractional number of primordial star-forming clouds that formed 
the eventual galaxy. However the associated population of IMBHs 
could amount to as much as ~ 10% of the current epoch stellar 
mass, if accretion onto these IMBHs also provides a possible ex- 
planation of the NIR background excess (Madau and Silk 2005). 
If such IMBH-induced jets carry flux, produced by a circum-black 
hole accretion disk dynamo, or drive jets that generate magnetic 
fields vis the proton Weibel instability in the shock interaction zone, 
then one might easily imagine seeding the clouds with a magnetic 
field that more than sufficed to allow MRI to operate. 

In summary, we have outlined two alternative pathways which 
dispense with the need for Population III to consist exclusively of 
very massive stars and for the transformation from first stars into 
Population II to be determined exclusively by the rise in gas and/or 
dust phase metallicity. MRI-dynamo generated magnetic flux in the 
protostellar accretion disk is the key to these possible scenarios. 
The nature of the seed field is the major difference. One option is 
that protogalactic Biermann battery-seeded MRI dynamos provide 
circumstellar disks with the necessary field strength to overcome 
high accretion rates characteristic of primordial enviroimients. This 
allows, and indeed the mechanism itself requires, stars of mass be- 
low 10 — IOOMq to form. In this case, there is no Population III. 
Another possibility is that the first objects (Population III is a mis- 
nomer since no nucleosynthetic tracers remain) were very massive 
( ^ IOOOM0), and form IMBHs which expel magnetic-flux-loaded 
jets that have been generated by dynamo activity in the circum- 
IMBH disk. The MRI prostellar dynamos are seeded by jet-induced 
instabiUties, and again, the first stars are expected to span a broad 
mass range. The IMF is a consequence of magnetic feedback. Of 
course, it may vary with epoch and/or environment, but the prin- 
cipal point we wish to emphasize is that it spans the conventional 
range of stellar masses. One could then make the transition from 
Primordial Star Formation to Current Epoch Star Formation while 
the mean overall metallicity remained extremely low. 
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